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ABSTRACT 

We present the analysis of the complete 3-year data set obtained with the Degree Angular Scale 
Interferometer (DASI) polarization experiment, operating from the Amundsen-Scott South Pole research 
station. Additional data obtained at the end of the 2002 Austral winter and throughout the 2003 season 
were added to the data from which the first detection of polarization of the cosmic microwave background 
radiation was reported. The analysis of the combined data supports, with increased statistical power, all 
of the conclusions drawn from the initial data set. In particular, the detection of -E-mode polarization is 
increased to 6.3a confidence level, TE cross-polarization is detected at 2.9a, and -B-mode polarization 
is consistent with zero, with an upper limit well below the level of the detected E-mode polarization. 
The results are in excellent agreement with the predictions of the cosmological model that has emerged 
from CMB temperature measurements. The analysis also demonstrates that contamination of the data 
by known sources of foreground emission is insignificant. 

Subject headings: cosmic microwave background — cosmology:observations — 
techniques:intcrferomctric tcchniques:polarimetric 



1. INTRODUCTION 

The use of the cosmic microwave background (CMB) as 
a precision probe of the composition and dynamics of the 
universe presumes the correctness of the theoretical frame- 
work by which we understand its origin and evolution. 
With the detections of CMB polarization anisotropy by 
DASI (Kovac et al. 2002, hereafter Paper V) and WMAP 
(Kogut et al. 2003) that framework itself has passed a 
critical test, as the presence of polarization is a generic re- 
quirement of the physics of recombination (for a recent re- 
view see Hu & Dodelson 2002). The intrinsic polarization 
of the CMB reflects the local radiation field anisotropy, 
specifically the quadrupole moment, at the surface of last 
scattering, 14 billion years ago. Within the context of the 
standard cosmological model, in which peaks in the CMB 
angular power spectrum are interpreted as the signature 
of acoustic oscillations seeded by primordial density fluc- 
tuations, current temperature measurements lead to spe- 
cific predictions for the shape and amplitude of the polar- 
ization power spectrum and the temperature-polarization 
cross power spectrum. Furthermore, it is a firm prediction 
that density fluctuations should create only _E-mode, i.e., 
curl-free, polarization patterns on the sky. 

In Paper V, we reported the high signal to noise (s/n) 
detection of i?-mode polarization in the CMB with DASI. 
Kogut et al. (2003) reported the detection of the TE corre- 
lation on larger angular scales by WMAP. The DASI and 
WMAP measurements underscore the power of combin- 
ing sub-orbital experiments, with their ability to achieve 
exquisite sensitivity over small areas of sky, with space ex- 
periments, whose sky coverage allows them to approach 
fundamental sample variance limits. These measurements 



provide fuel for a host of new experiments aimed at prob- 
ing the detailed polarization spectrum of the CMB, from 
the largest scales, which may bear the signature of in- 
flationary gravitational waves, to the smallest, where the 
CMB is expected to record the history of structure forma- 
tion. 

In Paper V, we presented the analysis of the initial DASI 
polarization data set obtained during 2001 and the first 
part of the 2002 season; modifications to the instrument 
to enable polarization measurements, as well as the ob- 
servations and calibration of the data, were described in 
Leitch et al. (2002a, hereafter Paper IV). In this paper we 
present the analysis of the 3-year DASI polarization data 
set, which comprises the initial data (271 days of observa- 
tions) and the additional data collected during the end of 
the 2002 season and over the entire 2003 season (191 days 
in all). The end of the 2003 season marked the conclusion 
of interferometric observations with DASI. 

2. OBSERVATIONS 

2.1. Instrument Review 

The use of interferometers, and of DASI in particular, to 
measure the angular power spectrum of CMB anisotropy 
is described extensively in Leitch et al. (2002b, hereafter 
Paper I) and in Paper IV. In short, an interferometer de- 
tects correlations between pairs of antennas, i.e., baselines, 
yielding measurements which are largely insensitive to un- 
corrected noise, and which map simply into Fourier com- 
ponents of the sky brightness distribution, longer baselines 
measuring higher spatial frequencies. As they operate 
natively in the Fourier plane, interferometers are ideally 
suited to measurement of the power spectrum. 
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DASI is a 13-element, co-planar interferometric array 
operating in ten 1-GHz bands spanning 26 — 36 GHz. The 
telescope sits atop an eleven meter tower located at the 
Amundsen-Scott South Pole research station. DASI was 
initially configured for measurements of the CMB temper- 
ature anisotropy on angular scales corresponding to mul- 
tipoles 140 < I < 900 (Paper I). These observations were 
completed during the 2000 Austral winter and the results 
reported in Halverson et al. (2002, hereafter Paper II) and 
Pryke et al. (2002, hereafter Paper III). As described in 
Paper IV, several changes were made for the polarization 
observations, including the installation of a large reflective 
radiation shield surrounding the telescope. The shield ge- 
ometry ensures that all lines of sight from the telescope 
in the direction of the station buildings, the horizon and 
ground, are reflected to the sky. To provide polarization 
sensitivity, the thirteen DASI receivers were each equipped 
with cooled, broadband achromatic polarizers which can 
be mechanically set to accept either left or right-circularly 
polarized light. Within each 1-hour interval of observa- 
tion, the polarizers were stepped between right (R) and 
left (L) positions to acquire data in all four Stokes states, 
which we refer to as co-polar (RR, LL) and cross-polar 
(RL, LR) data. 

A sun shield, designed to extend the polar observing sea- 
son to include periods just before sunset and after sunrise, 
was installed for the 2002 season. The shield was found to 
contaminate the shortest baselines by secondary reflection 
at times when the moon was above the horizon, and was 
removed for the final season. (Note however that the large 
ground shield remained in place during all three seasons.) 
Apart from the installation and removal of the sun shield, 
and replacement of miscellaneous hardware during routine 
maintenance of the telescope, the instrument was identical 
throughout the three seasons discussed here. 

2.2. Observing Strategy 

The polarization experiment consisted of deep observa- 
tions of two 3?4 FWHM fields, at Right Ascension 23^30 m 
and OO^O™, and Declination -55° (J2000) (see Figure 1). 
Details of the field selection are given in Paper IV. The 
fields were observed continuously at constant elevation 
(one of the unique advantages of observing from the South 
Pole), interspersed with daily observations of a calibra- 
tor source. A variety of cross-checks were built into the 
experiment design. Observations of the two fields were 
interleaved on 1-hour timescales to ensure that each was 
observed over the same azimuth range, allowing effective 
removal of any residual ground signals by differencing the 
data. The configuration of the array elements is three- 
fold symmetric and the entire array is fixed to a faceplate 
which can be rotated along the line of sight. On 24- hour 
timescales, the observations alternated between two differ- 
ent orientations of the faceplate, offset by 60°, in which the 
same Fourier components are measured, but by different 
pairs of antennas. 

Throughout the three seasons of observations, the ob- 
serving strategy, as described in Paper IV, remained iden- 
tical. 



2.3. Data Cuts 

A variety of cuts, discussed extensively in Paper IV, are 
applied to the data before input to the likelihood anal- 
yses; cut criteria and thresholds for the combined three- 
season data were the same as for the data presented in 
Paper IV and Paper V. During the 2002 season, when it 
was discovered that the short-baseline data were affected 
by reflections of the moon off the sun shield, a more strin- 
gent cut on the elevation of the moon was applied than 
during the 2001 season. During the 2003 season, for which 
the sun shield had been removed, cuts on the moon were 
relaxed to the 2001 levels. Early in the 2003 season, the 
polarizer assembly for one of the DASI receivers began to 
slip, and data associated with this receiver were therefore 
excised for the entire season, removing ~ 15% of the 2003 
data. 

Of the three-season data taken with the sun below the 
horizon, approximately 53% of the co-polar data, and 
61% of the cross-polar data passed the data cuts. Cut 
fractions were dominated by cuts on the elevation of the 
moon, and by cuts on the complex gains of multipliers in 
the correlator (see Paper IV). 

2.4. Calibration 

The instrumental response of each baseline of an inter- 
ferometer is the complex multiplication of the individual 
antenna responses. As a result, the polarization calibra- 
tion, up to an absolute phase offset, can be derived from 
observations of an unpolarized source (see Paper IV) . The 
remaining phase offset was determined by observing an 
unpolarized source through linearly polarizing wire grids 
attached to the front of each receiver. These offsets had 
previously been measured in 2001 August and 2002 Febru- 
ary. The wire-grid calibrations were performed again in 
2003 March and 2003 October, and results from all four 
epochs agree to within the measurement uncertainties. 

Leakage from one polarization state to the other will 
mix CMB power from T into E-modes and B-modes, and 
must be included in the analysis. Ideally these should be 
kept as small as possible, but even significant leakages can 
be accounted for if they are stable. The leakages can be 
separated into a component due to the polarizers (referred 
to as on- axis leakages in Paper IV), and a beam-dependent 
component which derives from the inherent properties of 
the optics (referred to as off- axis leakages in Paper IV). 
The latter were characterized in Paper IV, and are not ex- 
pected to change. The on-axis leakages for DASI were de- 
termined from observations of a bright unpolarized source 
in 2001 August, 2002 April, 2002 July and 2003 August. 
On-axis leakages of total intensity into the cross-polarized 
data were found to be less than 1% at all but the highest 
frequency bands (see Paper IV) and were consistent from 
year to year. 

3. RESULTS 

3.1. Noise Model and Consistency Tests 

For both the x 2 consistency tests and likelihood anal- 
yses that follow, the calibrated, field-differenced, and co- 
added visibility data are arranged in a data vector that is 
considered to be the sum of signal and noise components, 




Fig. 1. — Polarization field locations. DASI fields are marked by white circles, plotted over the Ka-band synchrotron map derived from 
the WMAP first-year data (Bennett et al. 2003). The color scale is linear and has been clipped at one thirtieth of the maximum to allow the 
structure at high Galactic latitudes to be seen. 



A = s + n. An accurate characterization of the instru- 
ment noise in the data vector is critical for estimating the 
faint polarization signal. We calculate various statistics 
from the raw data to test our assumptions about the noise 
model: namely, that the noise is Gaussian, that the noise 
variance is well approximated by the variance in the 8.4- 
s samples (the fundamental data rate) over 1-hr intervals 
before field differencing, and that the noise is uncorrelated 
between data vector elements. 

Details of the noise model tests and results for the initial 
data set are given in Paper V. We find similar results for 
the 3-year data. Three different types of short-timescale 
(1-hr) variance estimators are in excellent agreement with 
one another, and also with the variance derived from the 
scatter of 1-hr binned visibilities over the 3-year span of 
the data set. In Paper V, we scaled the noise by a factor 
of 1% (twice the maximum discrepancy among the various 
noise estimators described above) to test the robustness of 
our noise estimate and found no significant change in the 
likelihood results. 

We construct a data covariance matrix for each day of 
observation to test the assumption that the noise covari- 
ance matrix is diagonal, and to form a cut statistic based 
on correlated atmospheric noise. Data covariance matri- 
ces for all good-weather days are averaged, to test for noise 
correlations persistent throughout the data set. We find 
significant correlations in 0.017% of the off-diagonal terms. 
These are similar in number and magnitude to those pre- 
viously described in Paper V, where we found that adding 
the measured correlations to the noise model had a negli- 
gible effect on the final results. 

We construct a x 2 statistic on sum and difference data 
vectors for various splits and subsets of the data in which 
the same signal is expected, as a simple but powerful test 
of the noise model, data consistency, and presence of sig- 
nal in the data. The statistic is given by x 2 = A*Cl r A, 
where A = (Ax ± A2) /2 is the sum or difference data 
vector, and Cn = (Cjvi + CN2) /4 is the corresponding 
noise covariance matrix. We assess the significance of a x 2 
value by calculating the probability to exceed (PTE) that 
value in the cumulative distribution function. 

Following the procedure in Paper V, we split the data by 



epoch, azimuth range, faceplate position, and Stokes state 
(for co-polar observations) for each of the three years of 
data individually, for each year versus the other two, and 
for the full 3-year dataset. For each split, we also examine 
various subsets of the data vector, including co-polar and 
cross-polar data, data in distinct multipole ranges and the 
high expected s/n eigenmodes. In all, over 400 x 2 tests 
were performed on various splits/subsets of the data. 

The distribution of x 2 values for the differenced polar- 
ization data shows no evidence for residual signal in excess 
of the noise, giving a uniform distribution of PTE values 
with no apparent trends or outliers. In the 2003 tem- 
perature data, several subsets which isolate high expected 
s/n eigenmodes formed from the shortest baselines com- 
pared at different faceplate rotations yield high \ 2 values, 
the most significant of which gives a PTE = 0.0004. Be- 
cause this discrepancy appears in an isolated subset of data 
which contains an extremely high signal, it is suggestive 
of a small shift in instrumental response. Such a shift at 
the level suggested by these tests would not significantly 
affect the data. 

For the sum data vectors, a signal is detected with high 
significance in all of the co-polar splits and subsets, as 
well as in many of the high s/n subsets of the cross-polar 
data. In addition to the split data tests, we construct x 2 
statistics from the highest s/n eigenmodes of the unsplit 
data vectors, as a sensitive check for the presence of sig- 
nal, as discussed in Paper V. In the cross-polar data, the 
number of expected s/n > 1 polarization eigenmodes has 
increased from 34 in the initial data set to 60 for the 3- 
year data set (taking into account off-axis leakage terms in 
the noise covariance matrix Cat, see Paper V). For these 
eigenmodes, we find \ 2 = 141.3, which corresponds to a 
PTE of 1.66 x 10 -8 — a high-significance detection of a 
polarized signal. 

The results of both the x 2 consistency tests and the 
noise model tests demonstrate that the noise is well char- 
acterized and Gaussian, and that the data calibration is 
consistent over the three years of observations. The sum 
data vector x 2 tests on both the co-polar and cross-polar 
data indicate the detection of signal with high-significance, 
using a simple test that is independent of the more sophis- 
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ticatcd likelihood analysis results that follow. 

3.2. Systematic uncertainties 

In Paper V, we investigated the impact of possible 
sources of systematic uncertainties on the results of the 
likelihood analysis. In addition to the uncertainties in the 
noise model discussed above, we considered systematic un- 
certainties in the absolute cross-polar phase offsets, in the 
on-axis and off-axis leakages (see §2.4 for definitions), and 
in the pointing of the telescope. In each case, re-analyzing 
the data with the systematic errors set to the limits of their 
allowed ranges led to insignificant shifts in the results. 

The overall uncertainty of the DASI absolute calibra- 
tion, derived from measurements of external thermal loads, 
transferred to an astronomical calibrator, was estimated in 
Paper II to be 8% (ler), expressed as a fractional uncer- 
tainty on the Ci bandpowcr (4% in AT/T), and applies 
equally to the temperature and polarization data. 

3.3. Likelihood Results 

The same likelihood analysis formalism described in Pa- 
per V was used to analyze the DASI 3-year polarization 
data set. We place constraints on a given set of parame- 
ters by characterizing the shape of the likelihood function, 
i.e., the probability of our data vector A, given a set of 
parameters re which describe the signal, 

L(re) = P(A|re) 

oc dot (C (re))~ 1/2 cxp (~^A*C (re)" 1 A^j , (1) 

where C (re) is a model for the covariance matrix of the 
data vector. This matrix contains the expected contribu- 
tions of the instrument noise and the cosmological signal, 
as seen through the particular filter of the experiment. De- 
tails of the construction of C (re) for the DASI polarization 
data set are given in Paper V and in Kovac (2003). 

We report the maximum likelihood (ML) values of the 
parameters re, found using an iterated quadratic estima- 
tor technique (Bond et al. 1998). The shape of the like- 
lihood function in the vicinity of the ML peak is char- 
acterized using direct grid evaluation for those analyses 
with four or fewer parameters, and using a Markov chain 
algorithm for higher dimensional analyses (Christenscn 
et al. 2001). Confidence intervals for individual parame- 
ters are obtained by marginalizing (integrating) the likeli- 
hood function over the remaining parameters, and report- 
ing the equal-likelihood bounds which enclose 68% of the 
marginal likelihood distribution (the so-called highest pos- 
terior density (HPD) interval). These are given in the text 
in the format: ML (HPD-low to HPD-high). For param- 
eters which are intrinsically positive, we consider placing 
(physical) upper limits by marginalizing the likelihood dis- 
tribution as before, but excluding the unphysical negative 
values. We then test if the 95% integral point has a like- 
lihood smaller than that at zero; if it does, we report an 
upper limit rather than a confidence interval. We also in- 
clude in tabulated results the uncertainty estimates and 
parameter correlation coefficients obtained by evaluating 

1 A ACDM model with flat spatial curvature, 5% baryonic mat 
65 km s" 1 Mpc-\ {Cl b = 0.05, f2 cdm = 0.35, n A = 0.60, h = 0.65) 
ters changes the results of our T, E and TE amplitudes by less than 



the Fisher curvature matrix F at the peak of the likelihood 
function. 

Our goodness-of-fit statistic, as described in Paper V, 
is the logarithmic ratio of the maximum of the likelihood 
to its value for some hypothetical model Ho described by 
parameters kq, 

In the limit where the kml parameter uncertainties arc 
normally distributed, the statistic reduces to A = A% 2 /2. 
Where possible, we have tested the distribution of A using 
Monte Carlo simulations and found A% 2 /2 to be an excel- 
lent approximation. We therefore use the well-known \ 2 
cumulative distribution function to compute the probabil- 
ity that, assuming a given hypothetical model TCq to be 
true, we would obtain a value for A (Ho) that exceeds the 
observed value. This is quoted in the text and tabulated 
results as a PTE and, where appropriate, as a significance 
of detection, in units of a. 

For the 3-year data set, we perform all of the analy- 
ses reported in Paper V for the initial data set. In addi- 
tion, we perform two new analyses to constrain the level 
of contamination from synchrotron or radio point source 
emission. The separate likelihood analyses constrain dif- 
ferent aspects of the signal in the polarization data, in the 
temperature data, or in both analyzed together. The nu- 
merical results are given in Tables 1, 2 and 3, respectively. 
Closely following the format from Paper V, we now briefly 
discuss each analysis in turn. 

3.3.1. E/B Analysis 

In this analysis, two single-bandpower parameters are 
used to characterize the amplitudes of the E and B po- 
larization spectra. DASI has an instrumental response to 
E and B that is symmetric and nearly independent. Al- 
though the B spectrum is not expected to have the same 
shape as the E spectrum, we choose the same shape for 
both spectra to preserve this symmetry in the analysis. 

Following Paper V, we consider three a priori shapes to 
determine which is most favored by the data: the con- 
cordance 1 E spectrum shape, and two alternatives: a 
flat 1(1 + l)Ci cx constant spectrum, and a power law 
1(1 + 1)C; oc I 2 spectrum. For each of these three cases, 
the point at E = 0, B = corresponds to the same 
zero-polarization nopol model, so that the likelihood ratios 
A(nopol) may be compared directly to assess the relative 
likelihoods of the best-fit models. For the flat case, the 
ML flat bandpower values are E = 6.2/zAT 2 , B = 0.0/iK 2 , 
with K(nopol) = 10.58. For the I 2 power law case, the 
ML values are E = 4.2fiK 2 ,B = O.OfiK 2 (1(1 + l)Ci/2ir 
at I = 300), with A(nopol) = 15.41. For the concor- 
dance shape, the ML values are E = 0.73, B = 0.03 
in units of the concordance E spectrum amplitude, with 
A(nopol) = 20.10. Although each of these analyses shows 
strong evidence for an E polarization signal, the likelihood 
of the best fit model in the concordance case is a factor 
of 110 and 14,000 higher than those of the I 2 power law 
and flat cases, respectively. The data clearly prefer the 

', 35% dark matter, 60% dark energy, and a Hubble constant of 
id the normalization Cio = 700 fj,K 2 . Using the WMAPext paramc- 
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concordance shape, which wc therefore use for the E/B 
and other single bandpower analyses. 

Figure 2 illustrates the result of the E/B concordance- 
shape polarization analysis. The maximum likelihood 
value of E is 0.73 (0.53 to 0.95 at 68% confidence). For 
B, the result should clearly be regarded as an upper limit; 
95% of the B > likelihood (marginalized over E) lies 
below 0.25. 

The data are highly incompatible with the no polariza- 
tion hypothesis. The likelihood ratio A(nopol) = 20.10 
implies a probability that the data are consistent with the 
zero polarization hypothesis of PTE = 1.86 x 10 -9 , as- 
suming that the uncertainties in E and B are normally 
distributed, as discussed above. Marginalizing over B, we 
find A(E = 0) = 19.65, corresponding to detection of E- 
mode polarization at a PTE of 3.64 x 10~ 10 , or a signifi- 
cance of 6.27ct. 

The likelihood ratio for the concordance model gives 
A(E = 1, B = 0) = 0.68, for which the PTE is 0.51. We 
conclude that the data are consistent with the concordance 
model. 

3.3.2. E5/B5 

This analysis parameterizes the E and -B-mode spectra 
using five, piecewise-continuous, flat bandpowers for each. 
Results of this ten-parameter analysis are shown in Fig- 
ure 3. The Z-ranges of the five bands, provided in Table 1, 
are unchanged from the equivalent analysis presented in 
Paper V, with corresponding window functions similar to 
those previously reported. 

To test for consistency with the concordance model, 
we calculate the expectation value for that model for 
each of the five bands, yielding _E=(0.8,14,13,37,16) and 
B= (0,0,0,0,0) nK 2 . The ratio of the likelihood at this 
point in the ten dimensional parameter space to the max- 
imum likelihood gives A = 3.21, which for ten degrees 
of freedom results in a PTE of 0.78, indicating that our 
data are consistent with the expected polarization pa- 
rameterized in this way. The E5 / B5 results are also 
highly inconsistent with the zero polarization nopol hy- 
pothesis, with A = 19.47, which for 10 parameters yields 
PTE = 2.6 x 10~ 5 . This likelihood ratio is close to the 
value A(nopol) = 20.10 obtained for the best-fit E/B 
concordance-shaped model in the previous section, indi- 
cating that the two best-fit models, though obtained with 
different numbers of parameters, offer similarly likely de- 
scriptions of the polarization signal. 

3.3.3. E/0 E 

This analysis parameterizes the _E-mode polarization 
signal as above, as well as the frequency spectral index 
0e of this signal relative to the CMB (such that = 
corresponds to a 2.73 K Planck spectrum). As expected, 
the results for the .E-mode amplitude are very similar to 
those for the E/B analysis described above. The result 
for the spectral index is E = 1-39 (-0.21 to 2.92). The 
index is consistent with CMB, as A(0 E = 0) = 0.372 for 
a PTE = 0.39. It is inconsistent with synchrotron emis- 
sion, which is rejected at 2.65ct (A(0 = —2.9) = 3.52 for a 
PTE = 0.008). Although the constraint is not strong, it 
is nevertheless interesting in the context of assessing the 
level of possible contamination by non-thermal emission, 
as discussed in §4. 



3.3.4. E / sync 

Here and in the next section, we describe two new anal- 
yses using two-component models to constrain the level of 
foreground emission in the data set. In both of these anal- 
yses, the first component is the same CMB E power spec- 
trum as in the E/B concordance-shape polarization anal- 
ysis, with frequency spectral index = 0. For the E / sync 
analysis, the second component is defined to be a model 
of polarized synchrotron emission with E = B, flat power 
spectra (i.e., 1(1 + 1)C; oc constant) and a temperature 
spectral index = —2.9. The results of this analysis, given 
in Table 1 and shown in Figure 4, show that the amplitude 
of the CMB E component remains consistent with the con- 
cordance model, while the amplitude of the synchrotron 
component is consistent with zero. The 95% upper limit 
on the synchrotron component, E = B = O.QlpK 2 , is 
well below the level of the observed E signal, implying, as 
discussed in §4, that there is no significant contamination 
from foreground synchrotron emission in the data set. 

3.3.5. E/point 

A two-component model was also used to test the ex- 
tent to which the signal in our data could be due to point 
sources. For this E/point analysis, the second component 
is a model for polarized point sources, with E = B and 
1(1 + l)Ci oc I 2 power spectra, for which we have con- 
servatively assumed = —2. The results are given in 
Table 1 and shown in Figure 4. The amplitude of the 
CMB E component is again consistent with the concor- 
dance model, while the point source model component is 
consistent with zero. The 95% upper limit on the point 
source component, E = B = 0.98/iAT 2 at I = 300, is well 
below the observed signal level. We conclude that the pres- 
ence of point sources consistent with the data has at most 
a small effect on the polarization results (see also §4). 

3.3.6. Scalar/Tensor 

To constrain T/S directly from the polarization data, we 
perform a two-component likelihood analysis, using for one 
component the concordance E shape with B = for the 
scalars, and for the other component the predicted shapes 
of Et and Bt for the tensors. In principle, because the 
scalar B-mode spectrum is zero, this approach avoids the 
fundamental sample variance limitations arising from us- 
ing the temperature spectrum alone. However, the results 
of the E5/B5 analysis above indicate that we can place 
only upper limits on the E and B— mode polarization at 
the angular scales most relevant (I < 200) for the tensor 
spectra. It is therefore not surprising that our limits on 
T/S derived from the polarization spectra, as reported in 
Table 1, remain quite weak. 

3.4. Temperature Data Analyses and Results for the T 

Spectrum 

3.4.1. T/0t 

In this analysis, we use the bandpower shape of the con- 
cordance T spectrum, with the amplitude parameter ex- 
pressed in units relative to that spectrum. The spectral 
index is relative to the CMB, so that = corresponds 
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Table 1 

Results of Likelihood Analyses from Polarization Data 



68% interval 



analysis 


parameter 


/low 


— /high 


ML est. 


(F- 1 ) 1 / 2 error 


lower 


upper 


U.L.(95%) 


units 


E/B 


E 






0.73 


±0.21 


0.53 


0.95 




fraction o± concordance E 


B 






0.03 


±0.10 


-0.05 


0.12 


0.25 


fraction of concordance E 


E5/B5 


El 


28 


- 245 


0.06 


±0.8 


-0.58 


1.28 


3.6 


liK 2 




E2 


246 


- 420 


10.5 


±4.1 


6.3 


16.6 








E3 


421 


- 596 


3.6 


±4.7 


-2.1 


8.8 


15.2 


fj,K 2 




E4 


597 


- 772 


22.1 


±13.0 


9.8 


38.8 








E5 


773 ■ 


- 1050 


42.2 


±40.0 
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Fig. 2. — Results from the two-parameter shaped bandpower E/B polarization analysis. We assume the i?-mode power spectrum shape 
predicted for the concordance model, and the units of amplitude are relative to that model. The same shape is assumed for the B-mode 
spectrum, (right panel) The point shows the maximum likelihood value with the cross indicating Fisher matrix errors. Likelihood contours 
are placed at levels exp(— n 2 /2), n = 1, 2..., relative to the maximum, i.e., for a normal distribution, the extrema of these contours along either 
dimension would give the marginalized n-sigma interval, (left panels) The corresponding single-parameter likelihood distributions marginal- 
ized over the other parameter. Note the steep fall in likelihood toward low-power values; this likelihood shape (similar to a x 2 distribution) 
is typical for positive-definite parameters for which a limited number of high signal/noise modes are relevant. The grey lines enclose 68% of 
the total likelihood. The red line indicates the 95% confidence upper limit on .B-mode power. 



to a 2.73 K Planck spectrum. The amplitude of T has a 
maximum likelihood value of 1.23 (1.13 to 1.35), and the 
spectral index (3 T = 0.11 (-0.02 to 0.24). While the un- 
certainty in the temperature amplitude is dominated by 
sample variance, the uncertainty in the spectral index is 
limited only by the sensitivity and fractional bandwidth 
of DASI. Note that these results and those for T5 below 
have not been corrected for contributions from residual 
point sources (see Paper II). 

3.4.2. T5 

In the third panel of Figure 3, we present the results of 
an analysis using five, piecewise-continuous, flat bandpow- 
ers to characterize the temperature spectrum. Although 
these results are entirely dominated by the sample vari- 
ance in the differenced field, they are consistent with the 
far more sensitive temperature power spectrum described 
in Paper II. We include them here and in Paper V, pri- 
marily to emphasize that DASI makes measurements si- 
multaneously in all four Stokes parameters. 

3.5. Joint Analyses and Cross Spectra Results: TE, TB 

and EB 

3.5.1. T/E/TE 



In Figure 5, we present results of a three-parameter 
single-bandpower analysis of the T and E spectra, and the 
TE cross correlation spectrum, using bandpower shapes 
from the concordance model. The T and E constraints 
are quite similar to those from the E/B, E//3e and T / fir 
analyses described above, as expected. The additional re- 
sult here is the TE correlation, which has a maximum 
likelihood value of 1.08 (0.73 to 1.43). The data are in ex- 
cellent agreement with concordance model expectations. 

For the no X7?-correlation hypothesis, we find a likeli- 
hood ratio A( TE = 0) = 4.20 with a PTE of 0.0038; the no 
correlation hypothesis is rejected at 2.90c, a considerable 
improvement over the TE results presented in Paper V. 
Note that under the TE = hypothesis, negative and 
positive TE are equally likely, so the probability of falsely 
detecting a positive TE correlation at this level is PTE = 
0.0019, or 3.110-. 

3.5.2. T/E/TE5 

This is a seven-parameter analysis using single shaped 
band powers for T and E, and five piecewise-continuous, 
flat bandpowers for the TE cross correlation. The B-mode 
polarization is explicitly set to zero for this analysis. The 
T and E constraints are again similar to the values for 
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Fig. 3. — Results from several likelihood analyses. The top two panels show the ten-parameter E5/B5 polarization analysis, the third 
panel shows the T5 temperature analysis and the bottom panel shows the five TE bands from the T / E / TE5 joint analysis. All results are 
presented in flat bandpower units of 1(1 + l)C;/2-7r. The blue lines show the piecewise-continuous, flat-bandpower model for the maximum 
likelihood parameter values, with error bars indicating the 68% central region of the likelihood of each parameter, marginalizing over the 
other parameter values, analogous to the grey lines in the left panels of Figure 2. In each case, the solid green line is the concordance model 
used in Paper V and this paper, and the dashed dark-green line is a best-fit model for the WMAPext data from Spergel et al. (2003). 



the other analyses where these parameters appear. The considerable improvement from the results presented in 
interesting result here is the TE spectrum, shown in the Paper V, and trace the specific predictions of the concor- 
bottom panel of Figure 3. The bandpower results show dance model. 
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Fig. 4. — Results of E / 'sync and E/point two-component polarization analyses. These two analyses, performed to constrain the level of 
foreground emission in the data set, both use for their first component the same CMB E power spectrum as in the E/B concordance-shape 
polarization analysis (Figure 2). For the E/sync analysis shown in the upper three panels, the second component is a model for polarized 
synchrotron emission with E = B, flat power spectra and temperature spectral index /3 = —2.9 (see §3.3.4). For the E/point analysis shown 
in the lower panels, the second component models polarized point sources, assuming E = B spectra rising as 1(1 + l)Ci oc I 2 and j3 = -2 (see 
§3.3.5). The CMB E component is given in units of amplitude relative to the concordance model, as before. For comparison, the amplitudes 
of the second components are plotted in units relative to E = B = 8.5uK 2 at I = 300, which is the amplitude of the concordance model E 
spectrum at that angular scale. (The E/sync plot shows zero likelihood in the region where unphysical negative amplitudes for the second 
component cause the covariance to be non positive-definite). It is evident that the level of foreground contamination in either of these models 
is constrained to be well below the level of the observed E signal. 



3.5.3. T/E/B/TE/TB/EB 

Our last analysis is a six shaped-bandpower analysis for 
the three individual spectra T, E and B, together with 
the three possible cross-correlation spectra TE, TB and 
EB. Although there is no evidence for any £?-mode signal, 
we include the B cross-spectra for completeness. Since 



the TB and EB spectra are expected from theory to be 
zero, we preserve the symmetry of the analysis between E 
and B by parameterizing them in terms of the TE and E 
shapes, respectively. The results for T, E, B and TE are 
again similar to those obtained in the other analyses. As 
expected under the concordance model, there is no detec- 
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Table 2 

Results of Likelihood Analyses from Temperature Data 
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Table 3 

Results of Likelihood Analyses From Joint Temperature-Polarization Data 
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Fig. 5. — Results from the 3-parameter shaped bandpower T/E/TE joint analysis. Spectral shapes are as predicted for the concordance 
model, and the units are relative to that model. The layout of the plot is analogous to Figure 2. The two dimensional distribution in the 
right panel is marginalized over the T dimension. 



tion of EB or TB. 

4. DISCUSSION 

The analysis of the 3-year DASI data set supports all of 
the results presented in Paper V, in particular the detec- 
tion of -E-mode CMB polarization consistent with predic- 
tions of the standard cosmological model, and inconsistent 
with any significant contamination from foreground emis- 
sion. The increased sensitivity of the combined data set 
permits a more precise characterization of the detected 
polarization signal, and increasingly sensitive tests for the 
presence of foreground contamination. 

As discussed in §3.1, simple consistency tests with no 
assumptions about the underlying model show a high con- 
fidence detection of polarized signal. When the polar- 
ized anisotropy is modeled by a shaped spectrum over 
the /-range spanned by DASI, the data strongly prefer 
the concordance model shape for the -E-mode spectrum 
(see §3.3.1). Assuming the concordance model shape, we 



detect E-mode polarization (at 6.3cr confidence), as well 
as the TE cross-correlation (2.9cr), both with significantly 
increased statistical power over the results in Paper V, 
and at levels consistent with theoretical expectations. As 
expected for CMB anisotropies, the level of the B-mode 
polarization is consistent with zero, with an upper limit 
significantly below the E-mode detection. 

The five piecewise-continuous bandpower analysis pre- 
sented in §3.3.2 further demonstrates that the detailed 
shape of the E-mode anisotropy is consistent with the 
concordance predictions. The _E-mode bandpowers place 
a stringent upper limit in the lowest /-range band, with 
maximum likelihood values in the remaining four bands 
that are consistent with the damped concordance model, 
but inconsistent with a simple 1(1 + 1)C; cx I 2 power law. 
As shown in Figure 3 (§3.5.2), the TE angular power spec- 
trum provides a highly specific test of the spectral shape, 
as the predicted TE correlation crosses zero five times over 
the angular range covered by the DASI measurements. 
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Of the known foregrounds at 30 GHz, diffuse Galac- 
tic synchrotron emission and radio point sources are of 
the greatest concern. (As discussed in Paper V, neither 
free-free nor dust emission, of cither the thermal or spin- 
ning variety, is expected to contribute any significant po- 
larized signal to the DASI data.) It was for this reason 
that the two DASI fields were carefully selected from re- 
gions of exceptionally low emission in the Haslam maps 
(Haslam et al. 1981), confirmed at higher frequency by 
the WMAP synchrotron maps to be clean regions of sky 
(Bennett et al. 2003) (see Figure 1). The two fields were 
also selected from those with no detectable point sources 
from the original 32 DASI temperature anisotropy fields. 

While no published maps exist of the polarized syn- 
chrotron emission toward our fields, the DASI data them- 
selves set stringent constraints on the allowed level of dif- 
fuse synchrotron. As discussed in §3.3.4, a two-component 
model, designed specifically to probe for synchrotron emis- 
sion, finds that the amplitude of a synchrotron compo- 
nent with [3 — —2.9 is consistent with zero, with an up- 
per limit well below the level of the CMB i?-mode signal. 
The i?-mode spectral index constraint, improved from Pa- 
per V, indicates that the polarized signal is consistent with 
CMB; a signal with the expected synchrotron spectral in- 
dex (3 = —2.9 is rejected at the 2.7a level. The detection 
of the TE correlation is strong evidence that the polar- 
ized signal shares a common origin with the total intensity 
signal, a signal which has been robustly demonstrated in 
Paper II, Paper V and here to be consistent with CMB in 
both its frequency spectral index (/3 = 0.11 (—0.02 to 0.24) 
from §3.4.1 above) and angular power spectrum. More- 
over, the pronounced asymmetry of the _E-mode and B- 
mode levels allowed by DASI is in sharp contrast to the 
symmetric predictions for sources of foreground emission 
(Zaldarriaga 2001). 

Since no catalog of faint polarized sources exists over our 
observing region, and it is not possible to identify from to- 
tal intensity which point sources will have the strongest 
polarized flux, it is impracticable to project sources out 
of the polarization analyses. Nevertheless we can estimate 
the level of residual contamination in the data. A new, 
two-component analysis performed on the 3-year data set 
finds that the amplitude of a point-source component, with 
a conservative spectral index (3 = —2, is consistent with 
zero, with an upper limit well below the level of the CMB 
-E-mode signal (see §3.3.5). 

We analyzed maps of the differenced fields directly for 
evidence of individual sources appearing above the noise 
floor. The pixel value distributions are consistent with in- 
strument noise plus a Gaussian signal, and show no sign of 
excess tails due to discrete sources. The sensitivity of our 
maps is such that point sources with beam-attenuated po- 
larized flux greater than 15 mJy would have been clearly 
detected. 

In Paper V, we also presented a Monte Carlo analysis of 
the expected point-source contribution, based on available 
knowledge of the source populations. The number-flux 
counts used in the simulation (determined from the 2000 
DASI temperature data (see Paper V)) are consistent with 
the WMAP results on dN/dS (Bennett et al. 2003), and 



while there is a dearth of data concerning the distribu- 
tion of source polarization fractions at 30 GHz, a recent 
18.5 GHz study (Ricci et al. 2004) reports a distribution 
similar to the one assumed in our calculations. These sim- 
ulations suggest a mean bias of the E parameter of 0.04 
with a standard deviation of 0.05. We therefore conclude 
that there is neither the expectation of, nor evidence in 
the data for, any significant contamination by polarized 
point sources. 

5. CONCLUSION 

DASI has detected E-mode CMB polarization with high 
confidence (6.3a with the addition of the new data), and 
has detected the TE correlation at the 2.9cr level. The 
shape and amplitude of the detected i?-mode angular 
power spectrum are in agreement with theoretical pre- 
dictions based on the CMB temperature measurements. 
These results lend strong support to the underlying the- 
oretical framework for the generation of CMB anisotropy 
and point toward a promising future for the field — a fu- 
ture which includes the detailed characterization of the E- 
mode spectrum, the detection on small scales of B-modes 
from gravitational lensing, and the tantalizing prospect of 
constraining inflationary models via the detection of large- 
scale _B-modes from primordial gravitational waves. 

The DASI polarization results demonstrate the beauty 
of using interferometric techniques to provide well- 
matched filters for E- and B-mode measurements, and 
the inherent power of interferometers to reject uncorre- 
cted noise and control systematics. The DASI results 
would continue to improve with increased integration time, 
however the sensitivity required for the detection of the B- 
mode signatures demands large arrays, beyond the current 
scaling limitations of existing correlators. While there are 
currently no plans to upgrade DASI as an interferometer, 
advances in technology may once again make interferom- 
etry an attractive choice for pursuing the extraordinary 
control of systematics required by the next generation of 
CMB polarization experiments. 
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